Scale-Setting Using the Extended Renormalization Group and the Principle of 
Maximum Conformality: the QCD Coupling Constant at Four Loops 
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A key problem in making precise perturbative QCD predictions is to set tlie proper renor- 
malization scale of the running coupling. The extended renormalization group equations, which 
express the invariance of the physical observables under both the renormalization scale- and 
scheme- parameter transformations, provide a convenient way for estimating the scale- and scheme- 
dependence of the physical process. In this paper, we present a solution for the scale-equation 
of the extended renormalization group equations at the four-loop level. Using the Principle of 
Maximum Conformality (PMC) / Brodsky-Lepage-Mackenzie (BLM) scale-setting method, all 
non-conformal {/3i}-terms in the perturbative expansion series can be summed into the running 
coupling, and the resulting scale-fixed predictions are independent of the renormalization scheme. 
The PMC/BLM scales can be fixed order- by-order. As a useful reference, we present a systematic 
and scheme-independent procedure for setting PMC/BLM scales up to NNLO. An explicit 
application for determining the scale setting of R^+^-iQ) up to four loops is presented. By using 
the world average a^'^(Mz) ~ 0.1184 ± 0.0007, we obtain the asymptotic scale for the 't Hooft 
scheme associated with the MS scheme, A- 



conventional MS scheme. At 



'tH 
MS 

213ts® MeV. 



245ljQ MeV, and the asymptotic scale for the 
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I. INTRODUCTION 

All physical predictions in QCD should in principle be 
invariant under any choice of renormalization scale and 
scheme. However at any finite order, the use of differ- 
ent scales and schemes may lead to different theoretical 
predictions. The optimal procedure for obtaining precise 
QCD predictions is to choose the renormalization scale 
so that the result is scheme independent at any fixed or- 
der in as- The result of a scale-setting strategy should 
satisfy several self-consistent conditions: the existence 
and uniqueness of the scale, reflexivity, symmetry and 
transitivity Moreover, perturbative QCD becomes 
an Abelian theory as Nc — ^ 0, so QCD scale-setting must 
agree with that of QED in this limit [J]. We shall show 
that the Brodsky-Lepage-Mackenzie method (BLM) Q 
and the Principle of Maximum Conformality (PMC) [4j 
provide a solution to this problem ^ . 

The main idea of PMC/BLM is that after proper pro- 
cedures, all non-conformal {/3i}-terms in the perturba- 
tive expansion are summed into the running coupling so 
that the remaining terms in the perturbative series are 
identical to that of a conformal theory, i.e., the corre- 
sponding theory with {/J^} = {0}. The QCD predictions 
from PMC/BLM are then independent of renormaliza- 
tion scheme. It has been found that PMC/BLM satisfies 
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^ PMC provides the principle underlying BLM scale setting, so if 
not specially stated, we usually treat them on equal footing. 



all self-consistent conditions [1]. After PMC/BLM scale 
setting, the divergent "renormalon" series (n!/3f a") does 
not appear in the conformal series; thus as in QED, the 
scale can be unambiguously set by PMC/BLM. 

One can use PMC/BLM to relate perturbative calcu- 
lable observables in QCD, i.e. to derive commensurate 
scale relations among different observables [1, @ . More- 
over, from the requirement of scheme- independence, one 
can determine the displacements among the PMC/BLM 
scales that are derived under different schemes or conven- 
tions. We shall show how to fix the PMC/BLM scales 
order-by-order. The method for setting the leading-order 
(LO) and the next-to-leading order (NLO) PMC/BLM 
scales has been suggested in the literature [S'-'S]. In 
view of the recent improvements on perturbative QCD 
loop-calculations and the need to improve the theoret- 
ical predictions to confront more accurate experimental 
data, it is important to provide a systematic and scheme- 
independent treatment of PMC/BLM up to next-to-next- 
to- leading order (NNLO). 

We shall utilize a generalization of the conventional 
renormalization group (RG) analysis - extended RG 
equations which express the invariance of physical ob- 
servables under both the renormalization scale- and 
scheme- parameter transformations 0, In this ap- 
proach, a universal coupling function which covers all 
possible choices of scale and scheme is introduced, whose 
corresponding perturbative series serves as an intermedi- 
ate device for identifying the scale- and scheme- param- 
eters. It can be treated as a transparent solution to the 
scale-scheme ambiguity problem. A useful advantage is 
that the scheme dependence can be reliably estimated 
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through the scheme equations. This approach also pro- 
vides a platform for a reliable scheme-error analysis and 
gives a precise definition for the asymptotic scale under 
a possible renormalization scheme R, i.e., the scale for 
the 't Hooft scheme associated with the i?-scheme 

We shall present a general solution for the extended 
RG equation and give relations between the universal 
coupling function and the conventional adopted coupling 
function. 

The remaining parts of the paper are organized as 
follows: in Sec. II, we give the extended RG equa- 
tions and provide their solution up to four loops. In 
Sec. Ill, we present a systematic procedure for setting the 
PMC/BLM scales up to NNLO. Discussions and an ex- 
plicit application are also presented in Sec. III. Sec. IV pro- 
vides a summary. 



II. EXTENDED RENORMALIZATION GROUP 
EQUATIONS 

Conventionally, the scale dependence of an ordinary 
coupling constant is controlled by the RG equation or 
the /3-function 



simpler canonical form 



din /i^ 



47r 
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(1) 

where R stands for an arbitrary renormalization scheme. 
Various terms in /3(^, Pi' , • • •, correspond to one- loop and 
two-loop • • • contributions respectively. In general, the 
{/Sf"} are scheme-dependent and depend on the quark 
mass (m/) through the variable m^//i^. According to the 
decoupling theorem, the quark with mass m/ ^ /i can be 
ignored, and we can often neglect -terms when m/ <C 
fi. Then, for every renormalization scale /i, we can divide 
the quarks into active ones with nif = and inactive ones 
that can be ignored. Within this framework, it is well- 
known that the first two coefficients I3qi are universal, 
i.e., = 11 - 2n//3 and /3f = 102 - 38n//3 for Uf 
active flavors. Hereafter, we simply write them as /3o 
and /3i. Under the MS'-scheme, {I3]f^}i>2 up to four 
loops can be found in the literature Q. 

It will be convenient to use the first two universal co- 
efficients /3q and /3i to rescale the coupling constant and 
the scale-parameter In^^ in Eq.([T]). That is, by rescaling 
the coupling constant as 



Wo ^ 



(2) 

where cf = /3f /Sg" 7/3^ for i = 2, 3, • • •. 

As an extension of the ordinary coupling constant, one 
can define a universal coupling constant a(T, {ci}) to in- 
clude the dependence on the scheme parameters {cj 
which satisfies the following extended RG equations 



and 



da 



V [1 



020. 



(3) 



The scale-equation ([3]) , similar to Eq. ([2]) , can be used to 
evolve the universal coupling function from one scale to 
another. By comparing Eq.Q with Eq.Q, there exists 
a value oi t — tj^ for which 



(5) 



This shows that any coupling constant a^(T) can be ex- 
pressed by a universal coupling constant a(r, {q}) under 
the proper correspondence. The scheme-equation (|4]) can 
be used to relate the coupling functions under different 
schemes by changing {c^}. It is noted that the universal 
coupling function has a particularly simple form when all 
the scheme parameters {ci} are set to zero, i.e., the cou- 
pling function can be written as a function of the scale in 
terms of the Lambert W function 10]. This special case 
with {ci} = {0} is usually called the 't Hooft scheme 
[ll[. In addition to simplifying the solution of the RG 
equations, the 't Hooft scheme also provides a precise 
definition for the asymptotic scale (A) of QCD as will be 
shown below ^. 

The evolution of the universal running coupling can be 
obtained by integrating Eq.([3]), which can be rewritten 
as 



o 2 

Pi M5 



a(r,{c,}) 



da 



a(ro,{c.}) ^(a, {c»})' 



(6) 



where tq = (/3q//3i) In/Xg with /iq stands for an initial 
scale. Up to 0{a^), it leads to 



L = C+-+lna+(c2 - 1) a+ 



a +U(a ), (7) 



and the scale parameter as 

Pi 

one can express the RG equation H]) into the following 



^ Recently, it has been found that the 't Hooft scheme fails to re- 
produce the factorized form of the MS-scheme generalization of 
the generalized Crewther relation [T2| . This shows that one can- 
not use it for studying some special theoretical features of gauge 
theories beyond the two-loop level. Additional references and 
detailed discussions of the complimentary approach by Kataev 
et al. may be found in Ref. [T^ 
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where L = (/3q//3i) ln(/i^/A^) and C is an arbitrary in- 
tegration constant. A stands for the asymptotic scale, 
which is scale-invariant and leads to the coupling con- 
stant without any reference to the specific initial scale 
Mo [13, [3 • The value of A can be associated with the 
typical hadron size, which is not predicted by the QCD 
theory but must be extracted from a measurement of 
strong coupling constant at a given reference scale or a 
QCD measure with mass dimensions such as the pion de- 
cay constant The values of C and A are correlated 
with each other. One can find a general relation between 
the asymptotic parameters under different schemes from 
Eq.(l2l), i.e. for their values under two different schemes 
R and S, we have 



exp 



2^S 



As = exp 



/3i . ■ 



A, 



(8) 



where As and A/? are asymptotic scales corresponding to 
the choice of the integration constants Cs and Cr respec- 
tively. 

The 't Hooft scheme is free of higher-order corrections, 
i.e., all higher-order coefficients {ci} = {0}, and then its 
coupling constant is given by the solution of 



-'tH 



1 
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'tH 



In 
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where i'*^ = ln(^VA'*^^) and the integration 

constant C has been absorbed into the asymptotic scale 
A *^ for convenience. The 't Hooft coupling constant 



presents a formal singularity at L 



tH 



0: i.e. a 



tH 



a(0, {0}) — 00. Inversely, it can provide a precise defini- 
tion for the asymptotic scale; i.e., the 't Hooft scale A 



which is defined to be the pole of the coupling function in 
the 't Hooft scheme, a'*^ = a{^3^/|3l\n{^l'^/A''"^),{0}). 
Since the absorbed integration constant C is arbitrary, 
the value of A *^ is not unique, and there are infinite 
number of 't Hooft schemes, differing only by the value of 
A However, under a specific renormalization scheme 
(_R-scheme), its asymptotic scale can be fixed to be the 't 
Hooft scale associated with the _R-scheme A^^-^ Sj , which 
enters into both a^(^) = a(/3g//?i ln(^VAj?^^): {cf }) 



and 



'tH 



(m) 



'tH2 



),{0}). Here the 



word "associated" means we are choosing the particu- 
lar 't Hooft scheme that shares the same 't Hooft scale 
with the i?-scheme. In practice, one can obtain a relation 
between A^^ and the asymptotic scale A/j for i?-scheme 
by setting Cs — (here S stands for the 't Hooft scheme) 
in the left-hand side of Eq.®, i.e. 



exp 



As a special case, by choosing Cjjg 
we obtain 



Ar. (10) 



ln/?2//3i [T1[T1, 



MS \j2 



'/3i/2/3o' 



Ams' 



(11) 



which agrees with the observation presented in Ref.[8j. 
The present definition of Ajjg associated with the choice 
of C-jgg — In /3q //3i is the conventional one, originally 
suggest ed in Refs.il3l.[l4|: there are other choices for Cjjg 
[15l4l7l |. which can be helpful in certain cases. 

Eq.Q may be solved iteratively, and the solution to 
the universal coupling constant can be expanded as a 
power series of 1/L; i.e., up to four loops. 



^ + ^ {C - hiL) + ^ [C^ + C + C2 - (2C 



L2 



L3 

3c2 - 2 ) - 



I-C3 



InL-f l)lnL - l] 
I 



\tlL+-\ \tlL 



InL 



( ^ ) (12) 



As a cross-check, one finds that the above solution agrees 
with Ref. lis} after proper parameter transformations and 
by identifying the integration constant C* used there to 

heC* = ^{C- In ^) . When setting {cj = {0} and 
C = 0, we recover the coupling constant under the 't 
Hooft scheme. 



III. BLM SCALE-SETTING UP TO NNLO 

Generally, perturbative QCD prediction for a physical 
observable p can be written as 

p = ro [a:(Q) + (Ai + A2nf)a:+\Q) 
+{B^+B27if + B^n})a-+-'{Q) 
+(Ci +C2nf + Csn} + Cinj)al'+^{Q) + ■ ■ -(JlS) 

where as{Q) — ag{Q)/'K and the overall tree-level param- 
eter To is scale-independent and is free of as{Q). Here nf 
stands for the quark fiavor number and n{> 1) stands 
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for the initial as-order at the tree leveL After proper 
scale-setting, all n/-terms in the perturbative expansion 
can be summed into the running coupling. Here, we shall 
concentrate on those processes in which all ri/-terms are 
associated with the {/3i}-terms. Note that in higher-order 
processes, there are n/-terms coming from the Feynman 
diagrams with the light-by-light quark loops which are 
irrelevant to the ultra-violet cutoff. Those n /-terms have 
no relation to {/3i}-terms and they should be identi- 
fied and kept separate from the BLM scale setting ^ . 

The BLM scales can be determined in a general 
scheme-independent way. The generalization of the BLM 
procedure to higher order assigns a different renormaliza- 
tion scale for each order in the perturbative series. We 
can shift the rcnormalization scale Q into effective ones 
until we fully absorb those higher-order terms with n/- 
dependence into the running coupling LO and NLO 
BLM scale setting have been done in the literature [1, Q ■ 
Because of recent improvements in perturbative QCD 
loop-calculations, it is important to provide a system- 
atic and scheme-independent treatment of BLM up to 
NNLO. The BLM scales can be fixed order by order. In 
the following, we show how to set the BLM scales for the 
observable p. 

More explicitly, the first step of the BLM method is to 
set the effective scale Q* at LO 



P 



ro 



a^(Q*) + Aia':+'iQ*) + {Bi + B2nf)aJl+\Q*) 



+ {Ci + Can/ + C3n^)ar (Q*) + 



(14) 



The second step is to set the effective scale Q** at NLO 



p = ra 



+(Ci+C2n/)ar'(Q**) + 



(15) 



and the final step is to set the effective scale Q*** at 
NNLO 



To 



a';(g*)-f Aia';+^(Q**) + Sia';+^(Q***) 



+Ciar^(Q** 



(16) 



The step-by-step coefficients are presented in the Ap- 
pendix. When performing the scale shifts Q — Q* , 
Q* Q** and Q** — )■ Q***, we eliminate the n/-terms 
associated with the {/3i}-terms completely. At the same 
time, we also have to modify the coefficients, since the 



3 Those n/-terms, coming from the light-quark loops connected 
to at least four photon/gluon lines, are of higher twist and are 
usually power suppressed by the hard scales, so they can be safely 
neglected in typical applications. 

* Another wa y to set the BLM scale up to NNLO can be found 
in Refs.[19l |20( |. where however a unified effective scale Q* is 
adopted for all orders. 



net changes to the coefficients are proportional to /3- 
functions. To set the effective scale for a""'"'^, one needs 
even higher order information and here, a sensible choice 
is Q***, since this is the rcnormalization scale after shift- 
ing the scales up to NNLO. Note that the effective scales 
should be a perturbative series of as so as to absorb all 
n/-dependent terms properly, and up to NNLO, three 
effective scales can be written as 



In — = In — + — In —as{Q) + ( /^o In — 



*2 



?2, 2Q0 



-/?lln^)a2(g)^.o(a3^ 



(17) 



In- 



In 



Qo 



**2 



*2 ' 



^ln^a.(Q*) + 0(a?) (18) 



)*H=2 



**2 



+ Oias) (19) 



are determined so as 



where the effective scales Q^'** 

to eliminate A2nf, B2nf and C2n/-terms completely, the 
parameters x and z are used to eliminate the B^nj and 

the Can/ terms respectively, and the parameter y is used 



to eliminate the C4nf-term. It is found that 



Q2 



In 



**2 



Q 



*2 



In 



Q***2 


0**2 



6A2 

n 

6B2 
{n + l)Ii 

in + 2)Bi 



(20) 
(21) 

(22) 



and 



[n- 



nA2 
l)(2n-f 



(23) 



1)^3 _ e„(„ ^ 1)^2^3 + 6n^a 



nAl 



3(n 



2)Bl 



6(n+ l)AiC3 



(n + 1)^1^2 



(25) 



The coefficients Ai, Bi, Ci and etc. are renormalization- 
scheme dependent, so different rcnormalization schemes 
lead to different BLM scales g*:**,***; however the fi- 
nal result for p should be scheme independent. Using 
this argument, one can use BLM scale-setting method 
to relate perturbatively calculable observables; i.e. to 
derive commensurate scale relations among different ob- 
servables [Ej. In fact, any perturbatively-calculable phys- 
ical observable can be used to define an effective cou- 
pling constant by incorporating the entire radiative cor- 
rection into its definition [2^. For example Re+e-{Q) = 



-Rg+g- (Q) 1-1- defines an effective coupling con- 

stant Off (Q), where -Rg+g- (Q) stands for the Born result. 
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Any effective coupling constant can be used as a refer- 
ence coupling constant to define the renormalization pro- 
cedure. More generally, each effective running coupling 
constant or renormalization scheme is a special case of 
the universal coupling function as shown by Eq.([5]). 

The NLO commensurate scale relations between differ- 
ent effective coupling constants can be found in Ref.[5]. 
Replacing the observable p by its corresponding effective 
coupling constant and changing to be another effec- 
tive coupling constant, starting from Eq. p^ and follow- 
ing the same procedures, one can naturally obtain the 
commensurate scale relations up to NNLO. Moreover, by 
using the relations between Q* ■.**■*** and Q, one can find 
the needed scale displacement among the effective scales 
which are derived under different schemes or conventions 
so as to ensure the scheme-independence of the observ- 
ables. For example, from the relation between Q* and Q, 
one can easily obtain the well-known one-loop relation for 
the coupling constant [3], a^{e-^/^Q'^) = a^^^~^{Q'^), 
where the scale displacement e~^/'^ between the MS 
scheme and the Gell-Mann-Low scheme [23j is a result 
of the convention that is chosen to define the minimal 
dimensional regularization scheme (l3| . 



A. The PMC and BLM correspondence principle 

A systematic procedure for setting PMC scale at LO 
has been suggested in Ref.Q- The main procedure is to 
distinguish the nonconformal terms from the conformal 
terms by the variation of the cross section with respect to 
In /iQ (/io stands for some initial scale of the process). At 
LO, there is only one type of {/3i}-term (i.e. /3o) and the 
nonconformal terms always have the form of /3o In /ip , so 
one can determine the nonconformal terms exactly. How- 
ever, at higher orders, the In /ig-terms usually appear in a 
power series as In /Xg, (3i In /ig, P^iln /Zg)^, etc.. So this 
method is no longer adaptable to deal with the higher- 
order corrections, because the derivative with respect to 
a single In/ig cannot distinguish all the emerged {Pi}- 
terms. Some alternative should be introduced. 

The purpose of the running coupling in any gauge 
theory is to sum up all the terms involving the {ft}- 
functions, conversely, one can find all the needed {/3i}- 
terms at any relevant order from the expansion of the 
running coupling in the form of Eq. (jA9[) . Using this fact 
and also the known relation between {Pi} and Uf, one 
can obtain the PMC scales from the BLM scale-setting 
method. We call this the PMC and BLM correspondence 
principle. Since {ft} (i > 2) are scheme-dependent, the 



PMC and BLM correspondence depends on the renormal- 
ization scheme beyond the two- loop level. It is noted that 
such an expansion is different from that of Refs.[20l. [2l|. 
where all {ftj-terms which may contribute at the same 
order have been introduced to deal with the Adler D- 
function. 

More explicitly, up to NNLO, the physical observable 
p can be expanded in {ftj-series as. 



ro[a^(Q)-KA; + ^2/3o)«r'(Q) 

+ + C^^P2 + ClP^Pi + C°/33)<+3(Q)l .(26) 



The results for PMC can be naturally obtained from the 
BLM scale setting through proper parameter correspon- 
dence, i.e. 



A. = -l^i 

Bi = + 1025^ + 121B§ 
2 



B2 

B3 



3 
9 ^ 



{19B^ + 22Bl) 



(27) 
(28) 
(29) 
(30) 

(31) 



Ci = (7° 



2857 



+ 1122Cg + 1331C^ (32) 



C2 = - — (5033C2"- 3732C^- 4356C^) (33) 



C3 = — (325C2"-l-456Cg + 792C; 
C4 = -^Q 



(34) 
(35) 



which are obtained with the help of Eqs. (|13l26[) and the 
four- loop {/3i}-terms under the MS scheme ^]. 



B. An application of PMC/BLM scale setting 

We present an application of PMC /BLM scale setting 
up to NNLO by dealing with the total hadronic cross 
section in e+e" annihilation, R^+^~{Q) = R{e^e^ — J> 
hadrons). The explicit expression for Re+e-{Q) up to 
a^-order under the M^-scheme can be found in Ref . |24| . 
One finds 



i?e+e-(Q) =3^< 



(a^'HQ)) + A (a^''{Q))+ B (a^'^{Q)f + C (a^^^(Q))' 



(36) 
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TABLE I: Coefficients for the perturbative expansion of R^+^-iQ) before and after BLM scale-setting. 





Uf = 3 


Uf =4 


Uf = 5 


A 


1.6401 


1.5249 


1.4097 


B 


-10.2840 


-11.6857 


-12.8047 


C 


-106.8960 


-92.9124 +2«:/15 


-80.0075 +K/33 


A 


0.0849 


0.0849 


0.0849 


B 


-23.2269 


-23.3923 


-23.2645 


C 


82.3440 


82.3440 +2k/15 


82.3440 -I-k/33 



where 

A = 1.9857- 0.1152n/, 

B = -6.63694- 1.20013n/ - 0.00518?!^ - 1.240??, 

C = -156.61 + 18.777i/- 0.7974n^-f 0.02 1 57i^ + K7y, 

where r] = e,^ / ^3 J2q j is the electric charge 
for the active flavors. The coefficient k is yet to be deter- 



mined, its value is small [T^ and its contribution 

will be further suppressed by the factor 77, so we have 
set its value to zero in the following numerical calcula- 
tion. The values of A, B and C for = 3, 4 and 5 are 
presented in Tab HI 



At the present order in as, the n/-terms which come 
from the light- by- light quark loops and are irrelevant to 
the ultra-violet cutoff do not emerge, so all n/-terms in 



the above equation should be fully absorbed into ag ■ Af- 
ter doing the BLM scale setting up to NNLO, we obtain 



■e+e" 



(Q) 



+ (af ^(Q*)) + A (of ^(Q**))' + B (af ^(Q***))^ + C (af ^(Q***))' 



(37) 



where all the coefficients and effective scales can be cal- 
culated with the help of the formulae listed in the Ap- 
pendix. The coefficients are presented in TAB HI slight 

differences for B and C with varying Uf are caused by 
the charge-dependent parameter 77. 

From the experimental value, rp+p-iSLGGeV) = 
Y[Re+e- (31.6Gey) = 1.0527 ± 0.0050 we obtain 

Ag = 412;^°?MeV (38) 

and 

A^ = 359tlllMeY. (39) 

With the help of the four- loop coupling constant (IT^ . 
we obtain af^(Afz) = 0.129to:o?o- This value is some- 
what larger than the present world average af^^{Mz) = 
0.1184 ± 0.0007 dil. However, it is consistent with the 



values obtained from e+e colliders, i.e. a^^^{Mz) = 
0.13 ± 0.005 ± 0.03 by the CLEO Collaboration ^ and 

a¥^{Mz) = 0.1224 ± 0.0039 from the jet shape analysis 
|3l| . One may observe that a smaller central value of the 

world average for af'^{Mz) also results from the mea- 
surements of T-decays, T-decays, the jet production in 
the deep-inelastic-scattering processes, and from heavy 
quarkonia based on the unquenched QCD lattice calcu- 
lations A larger A^g^ leads to a larger af^'^(Mz), 
and vice versa. If we set af^^{Mz) to the present 
world average, we obtain A^|„^=5 = 245l?o MeV and 
A^|„,=5 - MeV ^ 



^ Ref . [32(1 obtained a slightly different value of ^'Ms\nj=5 = 
215 ± 9MeV. However, it is obtained by taking a wrong sign 
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It is found that after PMC/BLM scale-setting, the per- 
turbative expansion of Re+e-{Q — 31.6Gey) becomes 
more convergent. The higher-order corrections are used 
to set the scales Q* , Q** and etc., respectively. In par- 
ticular, we find Q* = (0.757 ± 0.008) Q which leads to 

a¥^{Q*)/a¥^{Q) = 1.060 ± 0.004. 

As a final remark, one can estimate the error caused 
by K with the help of the scheme-dependent equation 
(|3]). Such an analysis has been done in Ref.Q ®. It is 
found that even if we set its value that leads to the k- 
term has a comparable magnitude with those without k 
at the fourth-order, we shall only achieve an additional 
2% scheme error in addition to the above experimental 
errors. 



IV. SUMMARY 

The extended renormalization group equations pro- 
vide a convenient way for estimating both the scale- and 
scheme- dependence of the QCD predictions for a physi- 
cal process. The scheme dependence of a process can be 
reliably estimated by the scheme-equations for the ex- 
tended renormalization group. 

In the present paper, we have presented a general solu- 
tion to the scale equation of the extended renormalization 
group equations at the four-loop level. This formalism 
provides a platform for a reliable error analysis and also 
provides a precise definition for the asymptotic scale un- 
der any renormalization i?-scheme, A^^, which is defined 
as the pole of the strong coupling constant in the 't Hooft 
scheme associated with i?-scheme. 

We have also given a systematic and renormalization 

of (/33/2/3o) in the four-loop terms, which should be negative 
rather than positive. 
® Note there is a typo in Eq.(48) of Ref. 8], which should be 

changed to, ao = a+/^l + ^cf^a'^^ ^ . 



scheme- independent method for setting the PMC/BLM 
scales up to NNLO. The PMC provides the principle un- 
derlying ELM scale setting; the two methods are equiv- 
alent to each other through the PMC and ELM corre- 
spondence principle. The scales can be set unambigu- 
ously by PMC/BLM, which allows us to set the renormal- 
ization scale at any required orders to obtain a scheme- 
independent result. Such a scheme-independence can be 
adopted to derive commensurate scale relations among 
different observables and to find the displacements among 
the effective PMC/ELM scales which are derived under 
different schemes or conventions. 

The elimination of the renormalization scale ambiguity 
and the scheme dependence using PMC/ELM will not 
only increase the precision of QCD tests, but it will also 
increase the sensitivity of the collider experiments to new 
physics beyond the Standard Model. 
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Appendix A: Coefficients for the BLM scale setting 
up to NNLO 

The step-by-step coefficients for the BLM scale setting 
up to NNLO, which are introduced in Eas. (|14ll5ll6p . are 
listed in the following 



Co 



t 33 , 33 , ^ ^ 33 „ 

Ai + —A2 , Bi=Bi + —B2 , Ci = Ci + —C2 



1 

4n 



= — 



1 

4n 



1089(n + + \hZnA2 + 66(n + 1)^1^2 + (4Si - 1089^3)? 



66(n + \)A\ + 1971^2 + 4(n + 1)^1^2 - 4n(B2 + 33^3 



— ^-^ - 40392(74^3 143748A2''(3 + hn + 2n^) + 8A2n^(8Ci + 35937C4 
b4:A2n'^ 

5049B3n) - 13464:A2^n{n^ - 3n - 7) + 72^1^2(1 + n) (34^271 - 242^371 + 
121^2^(3 + 2n)) + 3A2^n(2857n -I- 352Bi(2 -I- n) - 95832^3(3 -1- 2n)) 



1 

192^2712 



22392C4n3 - 52272A2''(3 -I- 5n + 2n^) - 2AA2ii'^(-8C2 



(Al) 
(A2) 

(A3) 



(A4) 



8 



C2 



6534C4 + 933^371) - 48^1^2(1 + n)(19^2n - 13253^ + 66^2^(3 + 2n)) + 
^2^n-(-5033n - 192Bi(2 + n) + 316852(2 + n) + 52272^3(8 + 5n)) + 

24A2^n(-1871 + n(-627 + 311n)) 



1 



576^2^2 



2736C4n3 + 4752A2''(3 + 5?! + 2r?) + 144A2n2(4C3 + I98C4 + 



19B3?^) - 912A2''K - 4n) + 288^1^2(1 + n){-2B3n + Ai^S + 2n)) 



-.42^i(-325n + 57652(2 + n) + 9504^3(5 + 3n)) 



1 



4(n+ 1)^1 
-1 

4(n+ l)Ii 



33(n + 2)B2(2Bi + 33B2) + (n + 1)(153B2 + 4Ci - 1089C3)Ai 
2{n + 2)B2{2B'i + 335^) + {n + l)(19i?2 - 4.{C2 + SSCsMi 



I 

In deriving the above formulae, the following equation is can be obtained from its value at the scale Q, 
implicitly adopted, i.e. the value of Og at any scale Q* 



(A5) 



(A6) 
(A7) 
(A8) 



as{Q*) - as{Q) - ^/3oln ( ^ ) aj{Q) + 1 



Q*2 



/3iln(S 



1 

43 



Q* 



Q* 



aliQ) + 
ai{Q) + OiaUQ)). 



(A9) 
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